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An operational arrival-time distribution is defined as the 
distribution of detection times of the first photons emitted by 
two level atoms in resonance with a perpendicular laser beam 
in a time of flight experiment. For ultracold Cesium atoms 
the simulations are in excellent agreement with the theoretical 
ideal time-of-arrival distribution of Kijowski. 
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In many experiments the observables are the instants 
when certain events occur, or the durations of processes. 
This is also the case in experiments where the scale of 
the physical system requires a quantum treatment. Nev- 
ertheless, and, in retrospect, rather surprisingly, consid- 
ering time as an observable has been a long standing 
taboo for theorists because of Pauli's argument against 
the existence of a self-adjoint time operator conjugate to 
a Hamiltonian bounded from below ||l| , a condition that 
applies to most systems of physical interest. 

It is however quite clear nowadays, even though not 
yet a widespread piece of knowledge, that "Pauli's the- 
orem" is not a major stumbling block to formalize time 
in quantum mechanics, first because actual observables 
are in general linked to positive operator valued measures 
(POVM) that do not require self-adjoint (first-moment) 
operators 0, and second because in fact the "theorem" 
itself does not hold unless due care is taken of the do- 
mains of the operators implied (see for instance the anal- 
ysis in [|). 

While getting rid of the hindrance of Pauli's theorem 
is an important step, there still remains the question of 
defining time observables and/or time distributions. One 
then immediately runs into the difficulty of dealing with 
the "occurrence of events" in quantum mechanics, since 
the standard formalism provides a smooth, continuous 
description of the state evolution, without any explicit 
abruptness such as "collapses" or "quantum jumps" asso- 
ciated with the events. After all, this may not be surpris- 
ing if the quantum equations and the state are statistical 
and refer to ensembles; even clasically, abrupt changes for 
individuals (such as a yes/no transition) may be averaged 
out by statistical ensembles. The question then is how to 
extract information for individuals from the formalism, 
and in particular on the instants when the events occur. 
The current technical ability to trap and manipulate sin- 



gle atoms has triggered theoretical approaches to deal 
with a number of interesting time observables. These 
techniques make compatible the continuous and the jump 
descriptions, for the ensemble and the individual respec- 
tively. For example, the statistics of the dark periods of 
a fluorescing system with a rapidly decaying state and a 
metastable state have been theoretically understood and 
reproduced 0] . Our aim is to apply a similar analysis to 
arrival-time experiments. 

Among different time observables, the time of arrival 
(TOA) is one of the simplest, and an excellent study 
case since "arrival times" are routinely measured in par- 
ticle, nuclear, atomic, and molecular physics laboratories. 
There has been much recent interest among theorists on 
the arrival time, as reported in an extensive recent review 
g . In several works the possibility of defining ideal TOA 
distributions has been denied [|6|J2], or claims have been 
made that the arrival time cannot be measured precisely 
g. However, these negative conclusions have not been 
shared by all. In particular, Kijowski's TOA distribution 
for free motion has been defined and characterized as an 
optimal ideal distribution in the sense of satisfying in a 
unique manner several classically motivated properties 
H. Nevertheless, there is a clear divorce between the 
daily routine of many laboratories where time-of-flight 
experiments are performed, and the theoretical studies 
on the time of arrival, which are based on the parti- 
cle's wave function without recourse to extra (appara- 
tus) degrees of freedom. A number of "toy models" have 
been proposed that include simple couplings between the 
particle's motion and other degrees of freedom acting as 
clocks or stopwatches, but they do not incorporate any 
irreversibility and are still far from realistic experimental 
conditions 0. Halliwell has provided so far the only ir- 
reversible model for a time-of-arrival measurement based 
on a two level detector where the initial excited level 
decays due to the presence of the particle. The model 
remains however rather abstract and no identification is 
made with any specific measuring system |lO[ . 

Advancing in this direction, we provide in this letter 
a theoretical model to fill the gap between theory and 
experiment by capturing essential ingredients of current 
atomic time-of-fiight experimental settings, where the ar- 
rivals are measured as detections of fluorescence photons 
pT[ . Our model is based on a number of approxima- 
tions and idealizations, which are either standard ones, 
or can be easily improved upon to make it even more 



accurate and realistic. In the experiment we model, the 
atoms are first prepared in a certain state, then fly a dis- 
tance freely, and finally cross a perpendicular laser beam, 
resonant with one of the atomic transitions. The laser in- 
tensity is here assumed to increase sharply to a constant 
value in the direction of atomic propagation to facilitate 
the comparison with existing ideal TO A theories, and 
in particular with Kijowski's distribution. Note however 
that this "step" profile may be easily modified, e.g. by 
using a more realistic Gaussian shape. The first fluo- 
rescence photon detected marks operationally the arrival 
time since it is assumed that all emitted photons are de- 
tected by a perfect counter. The preparation is such that 
a single atom density operator will adequately describe 
the dynamics. Even though it is not strictly necessary for 
the model setting, we shall deal in the simulations with 
ultracold atoms at velocities of the order of centimeters 
per second. 

The atom's evolution can be modelled by a one dimen- 
sional, two component, quantum master equation. The 
rotating wave and electric dipole approximations, and a 
semiclassical description of the atom-laser coupling are 
assumed; the internal structure of the atom is reduced to 
the two levels coupled by the laser, |1) and |2), with zero 
detuning. The interaction picture based on the internal 
state atomic Hamiltonian Hint-, 



if,„, = i^(|2)(2|-|l)(ll). 



(1) 



is also used since it lets us eliminate the frequency u and 
the time explicitly from the Hamiltonian. Finally, assum- 
ing a vacuum state background (this is plausible because 
of the high value of optical frequencies) , the master equa- 
tion takes the form 

^ = -^[H,p] + l{2a-pa+ a+a'p pa+a'}, (2) 

where 



H = //2m + -n{x){a+ + a"). 



(3) 






|1)(2|, 
|2)(1|, 



fl{x) is the position dependent Rabi frequency (0(a;)ilo 
in the simulations), and 7 is the inverse of the lifetime 
of the excited state due to the coupling with the bath of 
modes of the background vacuum-fleld. 

Note that the master equation in Halliwell's model is 
formally similar pO[ | . There are however some important 
differences. In our case the two levels correspond to the 
internal structure of the particle, whereas in Halliwell's 
model they are assigned to the detector. Our atom be- 
gins in the ground state, it is excited by the laser and 
decays emitting a photon, whereas Halliwell's detector 
starts in an excited state and decays due to the coupling 



with the particle, which is restricted to a half-line. Halli- 
well's Hamiltonian in ph has only diagonal kinetic terms, 
and lacks the non-diagonal Rabi-frequency terms of (g) 
which are responsible for the absorption and emission in- 
duced by the laser. This leads in Halliwell's model to a 
closed dynamical equation for the population of the ex- 
cited state which does not hold here. Finally, while our 
dissipative term acts for all x, since an excited atom may 
also decay out of the laser field after the excitation, Hal- 
liwell's decay term is restricted spatially and is in fact 
the interaction term. 

The solution of the master equation (0) may be car- 
ried out by the techniques variously known as "quan- 
tum trajectories" 12|, "Monte Car lo wav e function ap- 



proach" (l3|, or "quantum jumps" |14| , [l5| . These meth- 
ods do not only allow to solve the master equation; they 
also provide, associated with speciflc measurement proce- 
dures, a temporal history or quantum trajectory for the 
individuals where the time observables can be read di- 
rectly. For our two level system, repeated measurements 
of the spontaneously emitted photons at time intervals 



dt < 



"0 : 



7 (so that at most only one spontaneous 



photon is emitted in dt), lead, on a coarser time scale, 
to a Schrodinger equation for the two component state 
vector conditioned to no photon detection, |'0c), 

{x\^/j,) = \l){x\i'c.i) + \2){x\iJc,2)- (4) 

This Schrodinger equation is governed by an effective 
Hamiltonian, 



H,ff^H-thl\2><2\ 



(5) 



The imaginary term responsible for the decay leads to a 
probability of detection of a spontaneously emitted pho- 
ton in the interval dt as given by Pc — "fdtP2 , where P2 is 
the probability of the excited state (note the caveats on 
identifying emitted and detected photons in |lq]). To re- 
produce the master equation, the ensemble is separated 
at every time step into detected or undetected cases by 
random collapses to the ground state according to pc [^ . 
Thus, each atom's history is characterized by a "trajec- 
tory" with a record of fluorescent photon detections at 
certain time steps. Averaging over many trajectories the 
master equation is recovered, but at variance with the 
smooth description of the master equation, each trajec- 
tory contains a sequence of collapses comparable to the 
experimental record of time-resolved photon counts for a 
single atom. For current purposes, only the time of detec- 
tion of the first photon is relevant, so in our case the time 
dependent norm of the state vector in the Schrodinger 
equation with the effective Hamiltonian Heff provides 
directly the operational arrival-time distribution H(t) we 
were seeking, as identical to the "delay function" P,p^, 

U{t) = -dN/dt, (6) 



7V = 



(^c,j(i)|^cj(t)), 



(7) 



J = l,2 



without the need of a random "Monte Carlo" sampling. 
This simplifies greatly the numerical treatment. 

In the simulations we have used the split operator 
method to solve the two-component Schrodinger equa- 
tion with the effective Hamiltonian da). For simplicity a 
minimum-uncertainty-product wave packet for a Cesium 
atom is prepared out of the laser range, with average ve- 
locities currently achieved by laser cooling, and with its 
center 1.05 iim away from the origin. Note that it is pos- 
sible to simulate longer flight paths since the free motion 
of the Gaussian wave-packet is known analytically. A 
closer modelling of specific experimental conditions may 
require combining several minimum uncertainty packets, 
for example to reproduce the spatially broad states gener- 
ated in a magneto-optical trap, or the direct propagation 
of initially mixed states rather than pure states. The ve- 
locity standard deviation, Sv = 0.098 cm/s, is taken in 
agreement with actual experimental conditions ll^l- We 
set 7 = 27r X 5.3 MHz as corresponds to the Cesium reso- 
nance line for the transition 6^6*1/2 ^^ &^Py,/2 at 852 nm, 
with a lifetime of 30.0 ns. 
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FIG. 1. The continuous line corresponds to Kijowski's 
distribution for the initial wavefunction described in the text. 
The longer dashed line is the operational time of arrival dis- 
tribution for Q.0 — 0.0997. Shorter dashes are associated to 
Q.0 = 0.3727, and dots to Q.o = I.247. 

The figure shows Kijowski's TOA distribution for the 
free motion of the initial state, with an average velocity 
of 10 cm/s, and the operational (normalized) TOA distri- 
butions obtained for three different values of JIq (0.0997, 
0.3727, and 1.247). For the smallest value the opera- 
tional distribution is clearly delayed with respect to Ki- 
jowski's, whereas for the two larger values the agreement 
is excellent. A remarkable finding is that the agreement is 
not broken by the significant fraction of atoms (9%) even- 
tually rejected by the laser field barrier for the largest fio- 
Essentially the same pattern with respect to the three 
values of fio has been also found for larger (40 cm/s) and 
smaller (1 cm/s) incident velocities. In the later case, 
only 20% of the atoms provide a fluorescence photon for 



f^o = 1-247. 

This model opens up many opportunities for further 
insight into other time observables, such as dwell times, 
or tunnelling times. In the case of the arrival time, it will 
stimulate the currently faint interaction between exper- 
iment and theory of time observables. It may also help 
the experimentalists to select beforehand optimal Rabi 
frequencies and lifetimes. For the theorist, it is a sim- 
ple tool to compare results of different approaches with 
simulations of experiments. An important bonus of the 
time distribution based on the first fluorescence photons 
is that it can be equally applied for free motion or in 
the presence of potentials (other than those associated 
to the measurement itself), so that it provides an opera- 
tional first-arrival-time distribution for which ideal (ap- 
paratus independent) theoretical distributions are not as 
well studied as for the free motion case. 
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